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We construct the general extension of the four-dimensional Jackiw-Rossi-Dirac Hamiltonian that 
preserves the antilinear reflection symmetry between the positive and negative energy eigenstates. 
Among other systems, the resulting Hamiltonian describes the s-wave superconducting vortex at 
the surface of the topological insulator, at a finite chemical potential, and in the presence of both 
Zeeman and orbital couplings to the external magnetic field. Here we find that the bound zero-mode 
exists only when the Zeeman term is below a critical value. Other physical realizations pertaining 
to graphene are considered, and some novel zero-energy wave functions are analytically computed. 



I. INTRODUCTION 

Eigenstates with precisely zero energy which are lo- 
calized in the core of topological defects in some order 
parameter in low-dimensional electronic systems have 
continued to draw attention. Historically, they were 
first studied in the context of fractionalization of elec- 
tric charge [3-[3]- Nowadays, they are also thought to 
provide a possible route to fault tolerant quantum com- 
putation 0t3| ; as well as a mechanism for ordering of the 
vortex core, for example (3-Hoj. Their existence is often 
tied to index theorems, which in turn rely on the sym- 
metry between positive and negative eigenvalues of the 
Hamiltonian. Such a reflection symmetry of the energy 
spectrum is certainly a necessary condition for the zero- 
energy states to be robust under variations in the defect 
potential: obviously, at least when their number is odd, 
at least one state has to be pinned to zero in order for 
the spectrum to remain symmetric under the exchange of 
the sign of the energy. This raises the following question: 
what is the general extension of the Jackiw-Rossi-Dirac 
3 (JRD) Hamiltonian (defined below in Eq. (1)) which 
still possesses the symmetry in question? In this paper 
we provide the answer and discuss some physical realiza- 
tions of our result. 

We show that irrespectively of a representation, there 
are exactly two, one linear and one antilinear, opera- 
tors which anticommute with the four-dimensional JRD 
Hamiltonian. The additional terms that respect the spec- 
trum's reflection symmetry and which can be added to 
it then have to be odd under one of these two operators. 
An important example of a Hamiltonian that anticom- 
mutes with the linear operator is the JRD Hamiltonian 
in presence of the Abelian gauge field [ll| ■ Here we focus 
on those terms that, in contrast, anticommute with the 
antilinear operator. In the four-dimensional representa- 
tion there are precisely four such terms. Their physical 
meaning is maybe most transparent in the context of the 
superconducting vortex at the surface of the topological 
insulator [12]: 1) the chemical potential, 2) the Zeeman 
coupling of the electron spin to the external magnetic 
field, and 3) the two components of the electromagnetic 
vector potential. Some special cases of our general Hamil- 
tonian have already appeared in literature in different 



physical contexts (l2| - [T3 | . It is shown here that the most 
general Hamiltonian with the vortex of unit vorticity in- 
deed has a zero-mode under a certain condition, and we 
exhibit the analytical solution for it when the flux of the 
gauge-field is localized at the origin. In the context of 
topological insulators, for example, we find that the Ma- 
jorana zero-mode [l2J exists only for the Zeeman term 
below a critical value, at which it delocalizes. 

The paper is organized as follows. In the next section 
we derive the antilinear operator that anticomutes with 
the JRD Hamiltonian. In Sec. Ill we display the ex- 
tended JRD Hamiltonian consistent with the antilinear 
spectral symmetry, and discuss its physical realizations. 
In Sec. IV we solve analytically for the zero-mode of the 
spectrum. The important example of a vortex at the sur- 
face of the topological insulator in the magnetic field is 
discussed in Sec. V, and we summarize our findings in 
Sec. VI. Some simple, but for our purposes crucial al- 
gebraic facts about Dirac matrices are presented in the 
Appendix. 

II. ANTILINEAR SPECTRAL SYMMETRY 

The JRD Hamiltonian 3] has the general form 

Ho = a\P\ + a 2 p2 + nii(x)a 3 + m 2 (2)0:4, (1) 

where on are four Hermitian four-dimensional matri- 
ces that satisfy the Clifford algebra {a,-, ay} = 2Sij, 
and the two masses are TOi(a?) = 2m(r) cos n<f> and 
rn,2(x) = 2m(r) sinrw/), with (r, </>) as polar coordinates in 
the plane, and the integer n is the vorticity. We assume 
m(r — > oo) to be finite, but the function m(r) other- 
wise arbitrary, pi are two components of the momentum 
operator. We have also set the velocity and the Planck 
constant to unity for convenience. 

It is well-known that there is a unique constant matrix 
that anticommutes with all four matrices Oj, and conse- 
quently with Hq: j3 — 01020304, which we chose here to 
be Hermitian as well. The existence of this matrix im- 
plies that the spectrum is symmetric around zero, and /3 
plays a crucial role in the derivation of the index theorem 
for the JRD Hamiltonian There exists, however, a 
unique additional antilinear operator that anticommutes 
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with Hq. To see this, the reader may recall that Hq may 
be understood as the effective low-energy Hamiltonian for 
a non-uniform, but real tight-binding Hamiltonian on the 
graphene's honeycomb lattice Irrespectively of the 
representation one can thus always define an operation 
formally (and in the case of graphene, literally |17| ) anal- 
ogous to time-reversal symmetry: an antilinear operator 
that commutes with Hq. Such an effective time-reversal 
operator may then be written as T — UK, where U is 
unitary, and K stands for complex conjugation. The ex- 
plicit form of U, of course, depends on the representation 
of the Clifford algebra. Since all four-dimensional repre- 
sentations are equivalent [l8|, , on the other hand, one 
is at liberty to chose the most convenient one for the 
present purposes. Let us assume therefore a representa- 
tion in which a.\ and ai are real, and cx 3 and on are imag- 
inary. The existence of such a representation is proved in 
the Appendix. In this representation the unitary part of 
T has to anticommute with all on , and therefore there is 
a unique solution, U — ft. The antilinear operator that 
anticommutes with Hq in this representation is then sim- 
ply 

A = @T = K, (2) 
where the first equality is completely general, and the 
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second is specific to the representation. We have used 
the fact that 1 = 1, and that the matrix j3 in this rep- 
resentation must be real. 

Uniqueness of the operator A can also be seen to fol- 
low from Schur's lemma [19j: since the unitary part of 
A would have to commute with all four a-matrices in 
the chosen representation, the irreducibility of the four- 
dimensional representation guarantees that it is propor- 
tional to unit operator. 



III. GENERAL HAMILTONIAN 

It is easy then to identify all four-dimensional purely 
imaginary Hcrmitian matrices that as such anticommute 
with the operator A. The first class consists of those 
which are even under T and anticommute with /3: a 3 
and «4, already present in Hq, and ia-i/3 = ia^a-iOt^ and 
ia2p = —iu\Ct-sa±, which are not. The second class con- 
tains terms that are odd under T and which commute 
with j3: ia,\Oii and ia^on. The most general Hamilto- 
nian, which while not anticommuting with (3 any longer 
still anticommutes with the antilinear operator A, and as 
such has the reflection symmetry of the spectrum around 
zero, may thus be written as 



H = ctj(pi + Ai{x){ia^ai)) + mi{x)a 3 + m2(x)a^ + \x{x)iazO.^ + h{x)ia\0t2- 

I 



(3) 



i=l 



It may also be useful to display the result using the stan- 7- matrices: ol\ = «7o7i, for i = 1, 2, 3, (X4 = ijojs'. 
dard covariant notation in terms of the Dirac (Hcrmitian) 
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H = ^2 Holiipi + ^(^735) + ijo(rnx(x)j 3 + m 2 (x)j 5 ) + (fi(x) + h(x)j )^ 35 . 

I 



(4) 



The matrix 735 = 17375 that plays the prominent role 
in H is the generator of the rotations of the two masses 
present in Hq into each other. Since 735 has the eigenval- 
ues ±1, Ai{x) enters the Hamiltonian as the axial gauge 
field, although we will see that it may also represent the 
true electromagnetic field in some physical realizations, 
notably in the topological insulator. h{x) is the mass- 
term which preserves the chiral symmetry of the mass- 
less Dirac Hamiltonian, generated by {73,75,735}, but 
breaks the (effective) time-reversal symmetry. 

The Hamiltonian in Eq. (4) represents the most gen- 
eral extension of the JRD Hamiltonian that preserves the 



antilinear symmetry between the positive and negative 
parts of the energy spectrum, with the additional terms 
independent of momentum. Momentum-dependent ex- 
tensions are also physically relevant, and will be a subject 
of a separate publication. (20j In this paper we have not 
considered the extensions of the JRD Hamiltonian that 
would preserve the linear spectral symmetry provided by 
/3 either. As already mentioned, an important example of 
the latter is the addition of the Abelian gauge potential 
by minimal subtraction, pi — > pi — A4, which represents 
graphene in magnetic field, for example, and also yields 
zero-energy states 11 1 . 
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Special cases of the Hamiltonian H have already arisen 
in different physical contexts. When [i(x) = h(x) = 0, 
Jackiw and Pi [HI have found that the axial vector po- 
tential can in certain sense be factored out of the Hamil- 
tonian, so that its presence does not alter the number of 
zero-energy states, but only modifies their form. When 
Ai(x) = h(x) = 0, the Hamiltonian describes the vortex 
in the s-wave superconducting order parameter at the 
surface of the to polo gical insulator at a finite chemical 
potential [lU, [2ll |22|] , and the exciton condensate in the 
symmetrically biased graphene bilayer [l3j . It may also 
be understood as describing the vortex in the Neel order 
parameter at finite Zeeman coupling to magnetic field in 
graphene, where one needs one copy of H for each Dirac 
point [8] . Similarly, the Bogoliubov - de Gennes Hamil- 
tonian for the vortex in a general superconducting order 
parameter on graphene's honeycomb lattice with the or- 
bital effect of the magnetic field included requires two 
copies of H with h(x) = [lOj. 

IV. ZERO-MODE 

Let us now proceed to solve for the zero-energy state 
for = 2/i and h(x) = 2h constant, when an ana- 

lytic solution is possible. In graphene, for example, we 
may chose a\ = —03 ® a%, a 2 = I ® 02, 0:3 = &i ® cri, 
and 0:4 = a 2 ® o~ 1 (23|. The Dirac fermion is then 
\I/ T = (tii, v%, U2, V2), with Ui and Vi standing for the 
components of the wave-functions at the two triangular 
sublattices of the honeycomb lattice, near the two Dirac 
points. (True spin degree of freedom is suppressed.) The 
masses are chosen so to represent the vortex in the Kekule 
bond-density- wave [l6| . The term proportional to \i is in 
this example the "pseudo-chemical potential" , which has 
the opposite sign at the two Dirac points, whereas the one 
proportional to h represents the pattern of imaginary- 
valued hoppings between the sites belonging to the same 
sublattice [24]]. The vector potential may be understood 
as arisingfrom a ripple on the graphene's surface, for 
example [17] . The equations for the zero-energy state in 
this representation become 

(/1 + h)ui + (id s - A)vi + fhv 2 = 0, (5) 

- ((j, + h)u 2 + mvi - (id z + A)v 2 = 0, (6) 
(ji - h)v 1 + (id z - A)u x + fhu 2 = 0, (7) 

- (fi- h)v 2 + mux - (idz + A)u 2 = 0, (8) 

where m = m(r)e 4 *, d z = e~ l ^(d r — (i/r)d^), A = 
A\ + iA 2 , and the bar denotes a complex conjugation. 
Assuming an ansatz 

e i{i-*){-)\ k = (-f+i g ( r ) e -S> (r ' )dr \ (9) 



e^ k+1 u k = f(r)e-Io m ^ dr \ (10) 
with k = 1,2, Eqs. (5)-(8) reduce to only two: 

(ji + h)g(r) + (a r +A 4 ,(r))f(r)=0, (11) 

(/x - h)f(r) -{d r + -- A4r))g(r) = 0, (12) 
r 

where we have also assumed a spherically symmetric field 
gauge field A(x) = A<p(r)(j). 

A. Zero gauge-field 

Let us first consider the problem without the gauge 
field, A$(r) = 0. Combining the last two equations, 
one finds the standard Besscl differential equation for the 
function g(r): 

r 2 d 2 r g{r) + rd r g(r) + [(p 2 - h 2 )r 2 - l]g(r) = 0. (13) 

For \h\ < the solution normalizable at the origin is 
qualitatively similar to the solution for h = [U [221 : 

g{r) = CJM» 2 -h 2 ) 1/2 ), (14) 

/(r) = C(^±|)VV (r( M 2 -/ l 2 ) 1 / 2 ), (15) 

where J n (z) are the Bessel functions of the first kind, 
and C is the normalization constant. The only difference 
from the solution at h = is in the characteristic length 
scale of oscillations, which now became longer. As long 
as > \h\ the zero-energy state is exponentially local- 
ized far from the vortex, with the characteristic length 
scale ~ l/m(oo). At fj, = h two of the components van- 
ish identically, v% — v 2 = 0, while the other two become 
constant. The zero-mode in this limit becomes the same 
as the one of the original JRD Hamiltonian Q. At the 
opposite end when fj, = —h, however, the solution is dif- 
ferent, f(r) — C, g(r) — C^r. In either case, at \fi\ = \h\ 
the characteristic length scale for the oscillations under 
the overall exponential decay of the solution diverges. 

When < I ft, I, on the other hand, the solutions turn 
into the modified Bessel functions 

g(r) = Ch(r(h 2 -^) 1/2 ), (16) 

f{r) = -C{^-^) 1/2 I {r(h 2 - /i 2 ) 1/2 ), (17) 
which now at large radius grow exponentially, 

9<r) « r ( fe2 _ M 2 )1/2 - (is) 



4 



The zero-energy state therefore remains normalizable 
only if the following condition is met: 

m(oo) 2 + /i 2 > h 2 . (19) 

When \h\ is above the critical value the zero- mode is ex- 
ponentially large far from the vortex. In a finite system 
this would presumably correspond to the state becoming 
localized at the boundary. Right at the critical value the 
zero-mode is thus critical: it oscillates with the amplitude 
of oscillations decaying as a power-law ~ 1/ s/r. 

B. Finite gauge- field 

With the general gauge-held present it seems not to 
be possible any more to find the analytic solution for 
the zero-mode, in contrast to the case when fi = h = 
We can nevertheless solve analytically a some- 
what artificial but nevertheless instructive example of the 
field's strength as the delta-function at the origin, when 
A^ir) = $/2r. To determine the sign and the magnitude 
of the flux $ that would render the energy of the vortex 
configuration assumed here finite, note that the Hamil- 
tonian H in Eq. (4) is invariant under the local unitary 
transformation, 

I 



One of the two constants appearing in the solutions is to 
be fixed by the condition at the origin, which, as usual, 
must be additionally specified [1H EEl , and which would 
correspond to different short-distance regularizations of 
the magnetic flux. 



Spreading the flux over a finite region of a linear size 
~ A around the origin may be accomplished by defining 
= (l/2A 2 )r for r < A, and A^ = l/2r, for r > A, for 
example. Very near the origin then the differential equa- 
tion would reduce to Eq. (13) with \i 2 — > + (1/2A 2 ), 
which would still yield one regular solution. At r = A the 
continuity of the solution would place two constraints on 
three (one for r < A, and two for r > A) constants of 
integration. The zero-mode clearly still exists. When 
\h\ > \fi\, modified Bessel functions replace those of the 
first kind far from the origin, and the zero-mode again 
delocalizes when the condition in Eq. (19) ceases to hold. 



provided it is accompanied by the gauge transformation 

A k -> A k + d k 6(it), (21) 

and the rotation of the complex mass, 

m^me 2ie(s \ (22) 

The last two transformations imply that the gauge- 
invariant coupling of the mass and the vector potential 
has the form 

\(d k -2iA k )m\ 2 , (23) 

so that the vortex configuration in the mass m has the 
finite energy only if A^ — l/2r at large distances. We 
thus take the flux to be $ = +1/2. Equations (11) and 
(12) then become quite similar: 

Oi + %(r) + (0 r + ^)/(r)=Q, (24) 



(/i - h)f(r) - & + -)g(r) = 0. (25) 

For | > \h\ the solutions are still the Bessel functions 
of the first kind, but now of the order ±1/2: 



(26) 



(27) 

I 

V. DISCUSSION 

The graphene representation we used is such that a.i 
are real for i = 1,3 and imaginary for i — 2,4. The 
time-reversal operator is then T = ia\a^K , and therefore 
A = ia2a^K = o\ ® a-sK. The zero-energy mode we 
found is therefore an eigenstate of the operator A with 
the eigenvalue +1. It is easy to check that the other 
eigenstate of A with the eigenvalue —1 has the opposite 
sign of the exponential in Eqs. (9) and (10), and therefore 
is not normalizable. 

The antilinear operator that anticommutes with the 
special case of our Hamiltonian when h = was already 
recognized in the second of the Ref. 14, in a specific rep- 
resentation similar to the one for graphene. The readers 
familiar with the BCS theory of superconductivity may 
also recognize it as being closely related to the ubiquitous 
symmetry of the BCS-type Hamiltonians, which origi- 



g(r) = C x J 1/2 (r(/i 2 - h 2 ) 1 ' 2 ) + C 2 J_ 1/2 (r(/i 2 - h 2 ) 1 ' 2 ), 
fir) = (^) 1/2 [Cr J 1/2 (?V - ^ 2 ) 1 / 2 ) _ C 2 J„ 1/2 (r( M 2 - h 2 f' 2 ) 
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nates in the Bogoliubov-Valatin doubling of degrees of 
freedom, characteristic for the BCS problem. We saw 
here, however, that the antilinear spectral symmetry is 
in fact a general property of certain wide class of Dirac 
Hamiltonians, which derives from the properties of the 
representations of Clifford algebra under complex conju- 
gation. As such it is also present in the realizations of 
the Dirac Hamiltonian relevant to graphene, where the 
masses may equally represent insulating order parame- 
ters [H [23 . 

At the surface of the topological insulator, we can con- 
sider the BdG Hamiltonian in presence of the vortex in 
the s-wave superconducting order parameter [l2l . l28j , by 
constructing the Dirac fermion as , 5 T = (c-j-, cr,, c|, ct). If 
the single particle Hamiltonian is h = p\o~\ + P2&2 ~ Mi 
the BdG Hamiltonian at fi — assumes the form in Eq. 
(1) but with the a-matrices as u\ = I®o~i, S2 = 03(g) 02, 
ct3 = eri<8>(72, «4 = <72<8>er2- Although this representation 
appears rather different from the one we used, the two 
are, of course, equivalent: cti = U^onU, for i = 1,2,3,4, 
with the unitary operator as 



VI. CONCLUSION 

In conclusion, we have determined the general exten- 
sion of the four-dimensional Jackiw-Rossi-Dirac Hamilto- 
nian that retains the antilinear reflection symmetry of the 
spectrum, and solved for its zero-energy state in several 
examples. A particularly relevant physical realization of 
the most general Hamiltonian of this kind is provided by 
the superconducting vortex at the surface of a topolog- 
ical insulator with the vortex, and/or an external mag- 
netic field fully included. The Majorana fermion inside 
the vortex core is found to exist only when the Zeeman 
coupling of electrons to the magnetic field is sufficiently 
small. 
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U = e l ^® a ^[I ® e l ^][<J3 ® da]- (28) 



Since the matrix —03 ® I = ia^a^ is the particle number 
operator, both the chemical potential and the electro- 
magnetic field enter the BdG Hamiltonian by multiplying 
it, just like in Eq. (3). Finally, as the generator of spin 
rotations around the direction perpendicular to the plane 
of the system is 03 <g> 03 — —ia\a,i-, the Zeeman coupling 
of the electrons to the magnetic field of the vortex, en- 
ters the BdG Hamiltonian precisely as the last term in 
Eq. (3). 

The superconducting vortex at the surface of the topo- 
logical insulator provides therefore a physical realization 
of the most general Hamiltonian with the antilinear re- 
flection symmetry of its spectrum. As the magnetic 
field always accompanies the superconducting vortex, we 
conclude that the Majorana fermion in the vortex will 
survive only if the Zeeman coupling of the electrons in 
the topological insulator to the magnetic field is small 
enough. In this respect we may note that our assumption 
of a constant Zeeman term is a reasonable approximation 
in the strong type-II limit, in which the magnetic field de- 
cays over the length scale of penetration depth A, whereas 
the zero-energy state decays over the much shorter super- 
conducting coherence length, £ ~ l/m(oo). According to 
Eq. (19), however, even for an overly strong Zeeman cou- 
pling the Majorana fermion can always be produced in- 
side the vortex by simply increasing the chemical poten- 
tial. Finally, one can imagine placing the whole system 
in an additional uniform magnetic field, which can then 
be used to manipulate the localization of the Majorana 
zero-mode. 



VIII. APPENDIX 

If one demands that five four-dimensional anticommut- 
ing Hermitian matrices, such as cti i — 1, ..4 and /3, are 
all either real or imaginary, it is easy to see that in any 
representation precisely three of these will be real and 
two imaginary. 

A direct way to show this is to construct the matrices 
cti and (3 out of standard Pauli matrices. The first three 
one may chose to be 

<7fe ® <7l,(7fc ® cr 2 , erfc ® 0-3 • (29) 

The remaining two will then be 

a n ® I, a m ® I (30) 

where n 7^ m 7^ k. Choosing k — 1 or k — 3 then makes 
the second in the first group and one in the second group 
imaginary, and the rest real. Choosing k — 2, on the 
other hand, makes the first and the third in the first 
group imaginary, and the rest real. Since all sets of five 
four-dimensional anticommuting Hermitian matrices of 
definite symmetry under complex conjugation are either 
exactly like in this example, or with the the two factor 
spaces interchanged, this proves our assertion. 

Another way to prove it would be to notice first that it 
would be impossible to have three of the matrices imag- 
inary and two real, since that would contradict the fact 
that any of the five matrices is a product of the remain- 
ing four. So a priori the only other options would be 
to have all five matrices real, or four imaginary and one 
real. Both cases would imply that there exists a four- 
dimensional representation of the Clifford algebra of four 
elements that all square to +1 or all to -1, which is purely 
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real, since in the latter case we could pick the four imag- 
inary matrices and multiply them by the imaginary unit. 
But that would, on the other hand, be in contradiction 
with the result that the smallest real representations of 
Clifford algebras C(4, 0) (four anticommuting elements 
each squaring to +1) and C(0, 4) (four anticommuting el- 
ements each squaring to -1) is actually eight-dimensional 



[29[ . The only Clifford algebras of four elements that ac- 
tually possess a real four-dimensional representation are 
C(3, 1) and C(2, 2), in accord with our result. 

In our derivation of the antilinear symmetry of the 
JRD Hamiltonian wc could indeed therefore pick two real 
matrices for those multiplying the momentum operator, 
and the two imaginary for the masses. 
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